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ABSTRACT 

We model the reflected spectrum expected from localized magnetic flares above 
an ionised accretion disk. We concentrate on the case of very luminous magnetic flares 
above a standard accretion disk extending down to the last stable orbit, and use a 
simple parameterisation to allow for an X-ray driven wind. Full disk spectra includ- 
ing relativistic smearing are calculated. When fit with the constant density reflection 
models, these spectra give both a low reflected fraction and a small line width as seen 
in the hard spectra from Galactic Black Hole Binaries and Active Galactic Nuclei. We 
fit our calculated spectra to real data from the low/hard state of Nova Muscae and 
Cyg X-1 and show that these models give comparable to those obtained from the 
constant density refiection models which implied a truncated disk. This explicitally 
demonstrates that the data are consistent with either magnetic fiares above an ionized 
disk extending down to the last stable orbit around a black hole, or with non-ionized, 
truncated disks. 
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1 INTRODUCTION 

X-ray reflection is in principle a sensitive diagnostic of the 
X-ray source geometry. The amount of reflection depends 
on the solid angle subtended by the optically thick material, 
(and also on the ionisation and elemental abundances in the 
reflector) while the amount of relativistic smearing of the 
reflected spectral features (especially the iron Kq line and 
edge) shows how far this material extends down into the 
gravitational potential (see the review by Fabian et al. 2000 
and references therein). This is important as there is much 
current debate about the structure of the X-ray source. One 
of the possibilities discussed in the literature is magnetic re- 
connection above the disk (e.g. Galeev, Rosner & Vaiana 
1979, Haardt, Maraschi & Ghisellini 1994) while another is 
the advection dominated accretion flows (e.g. Narayan & 
Yi 1995). These result in very different source geometries 
- in the advective flows the inner disk is truncated, being 
replaced by an X-ray hot plasma, while the magnetic recon- 
nection models have a disk which extends down to the last 
stable orbit. 

The reflected spectrum is generally calculated assuming 
a given density structure of the disk (generally constant) as 
a function of height, and with a constant ionisation state 
as a function of radius. Fits with these models to spectra 
from both Galactic Black Hole Candidates (GBHC) and 



some Active Galactic Nuclei (AGN) show overwhelmingly 
that the solid angle subtended by the disk is generally less 
than unity, and that the relativistic smearing is less than ex- 
pected for a disk which extends down to the last stable orbit 
(Zycki, Done & Smith 1997; 1998; 1999; Zdziarski, Lubihski 
& Smith 1999; Chiang et al., 2000; Done, Madejski & Zycki 
2000; Miller et al., 2001). They also show that the soUd angle 
and amount of smearing are correlated with spectral index 
in the sense that harder spectra show a smaller amount of 
reflection (f2/27r << 1) and narrower line emission (Zycki 
et al., 1999; Zdziarski et al., 1999; Gilfanov, Churazov & 
Revnivstev 1999; 2000; Lubihski & Zdziarski 2001). Ionisa- 
tion of the reflected spectrum is only strongly present for 
steep spectrum sources. 

This rules out the simplest version of the magnetic flare 
model unless the flares generate a mildly relativistic outflow 
away from the disk (Beloborodov 1999). However, such out- 
flows need not be required if the modelled reflected spectra 
are systematically underestimating the amount of reflection 
present in the data. For example, if the inner disk is so highly 
ionised that even iron is completely stripped, then reflection 
produces no spectral features (e.g. Fabian et al., 2000 and 
references therein). However, the data require a fairly sharp 
transition between the extreme ionisation and mainly neu- 
tral material, so as not to produce an intermediate region of 
high ionization, where iron is mainly H- and He-like, which 
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produces huge and unmistakable spectral features which are 
not present in the data (Done & Zycki 1999; Done et al., 
2000; Young et al., 2001). 

It has recently been realised that X-ray illumination 
changes the density structure of the disk. The X -ray heated 
material at the top can expand outwards, lowering its den- 
sity and increasing its ionisation. Calculating the self con- 
sistent density (and hence ionization) structure is espe- 
cially important as there is a thermal ionization instability 
which affects X-ray illuminated material in pressure bal- 
ance, which can lead to a hot, extremely ionised skin form- 
ing on top of the rest of the cooler, denser, mainly neutral 
disk material (Field 1965; Krolik, McKcc & Tarter 1981; 
Kallman & White 1989; Ko & Kallman 1994; Roiariska & 
Czerny 1996; Nayakshin, Kazanas & Kallman 2000, here- 
after NKK; Rozaiiska & Czerny 2000; Nayakshin & Kallman 
2001, hereafter NK; Ballantyne, Ross & Fabian 2001). 

The existance of the thermal instability can cause a sub- 
stantial underestimate of the amount of reflection present 
as measured by the fixed density reflection model spectra 
(NKK; NK). This is exphcially shown by Done & Nayak- 
shin (2001, hereafter DN) and Ballantyne et al., (2001) for 
reflection from a single radius in the disk. Here we extend 
the calculations of DN to build full disk reflection models for 
the case of magnetic flares. We give a simple parameteriza- 
tion of the one major uncertainty in these models, which is 
the strength of any sideways expansion of the X-ray heated 
material from underneath the flare (see §2). We first fit the 
calculated total spectra (incident spectrum plus reflection 
from a disk with self consistent density structure including 
the thermal instability) with the power law plus fixed den- 
sity reflection models in §3. We then fit the new models to 
actual data in §4 and show for the first time that the mag- 
netic flare models (in which the disk subtends a solid angle of 
27r) can fit the observed apparently weak reflection/narrow 
line, flat spectrum GBHC rather well. However, with these 
models, the data reqmre moderate outflow rates from the 
X ray heated region. Future work should explicitly calcu- 
late this sideways expansion of material, so as to remove the 
one remaining free parameter from the simplest magnetic 
flare model. 



2 REFLECTED SPECTRA 

The code of NKK gives the reflected spectrum from an X- 
ray illuminated slab whose vertical structure is given ini- 
tially by that of an accretion disk (Shakura & Sunyaev 
1973) at radius, r, and dimensionless mass accretion rate, 
m = riMc^ /IjEdd = L/LEdd (IjEdd is the Eddington accre- 
tion rate, L is the disk bolometric luminosity, M is the mass 
accretion rate and r] = 0.06 is the radiative efficiency of the 
standard disk in Newtonian limit), where a fraction f corona 
of the energy liberated by this infalling material is dissipated 
in a corona rather than in the optically thick disk (Svensson 
& Zdziarski 1994). This also defines the soft flux Fdisk at 
the top of the disk due to viscous dissipation underneath. 
There is also additional heating from the illuminating X ray 
flux _Fx with spectral shape Fx(E), which then changes the 
vertical structure since the heated upper layer can expand. 
The steady state density structure is calculated under the 
hydrostatic equilibrium assumption. 



The reflected spectrum from a full disk is an integral 
of reflected spectra from different radii. In practice we com- 
pute spectrum for only 6 radii due to CPU time limitations: 
r = 3.5, 4.9, 7, 14, 35 and 105 Schwarzchild radii for the 
same m and X-ray spectral shape, with a given radial de- 
pendance of the X-ray illuminating flux which depends on 
the X-ray source geometry. If the X-ray source(s) are mag- 
netic flares due to the MHD dynamo responsible for the vis- 
cosity (Balbus & Hawley 1991) then the X-ray flux might be 
expected to scale with the disk flux i.e. Fx oc Fdisk cc r^"^ so 
that Fx/ Fdisk is constant. The spectra from different radii 
can then be convolved with the relativistic transfer func- 
tions (assuming a Schwarzchild metric, but not including 
light bending; Fabian et al. 1989), which yields the requisite 
full disk relativistically broadened spectra. 

Given the four initial parameters (m, fcorona, Fx/Fdisk, 
and X-ray spectral shape), together with the assumption 
that Fx (X Fdisk, the reflected spectra from magnetic flares 
should be completely determined for a given inclination 
angle. However, the assumption of hydrostatic equilibrium 
is not entirely justified because a local outflow is formed 
(NKK, NK). Hydrostatic balance predicts that the ver- 
tical scale height for a skin of temperature T is Hs ~ 
{'ikTR^ /GMnipY^'^ . Compare this with the standard disk 
vertical pressure scale height, H, using equations of Svens- 
son & Zdziarski (1994): 

= 2.5 X 10-^Tr ( ^^^ 'o.o7°"°0 ^ ' 

where Tf is the skin's temperature in units of 10^ Kelvin, 
rh is the dimensionless accretion rate, r = R(? /2GM, and 
J(r) = 1 — \/3/r. At radius r = 6, for example, this ratio 

becomes 

(§)^6.5T.( -^(^-^7-y \ (2) 

For Cyg X-1, rh(l — /) is likely to be less than 0.01, and 
Ty ^ few for hard X-ray spectra. Now, magnetic flares are 
likely to have height not much larger than ~ H (Nayakshin 
& Kazanas 2001b) above the disk. Hence, the hydrostatic 
balance condition predicts that the skin can extend to height 
Hs ^ H, i.e., completely surround the X-ray source. This is 
inconsistent with the explicit assumption of the NKK code 
that X-rays are incident on the top of the skin at a given 
angle (equivalent to the source located at infinity from the 
disk). 

In addition, when the X-ray source is engulfed by the 

highly photo-ionized skin, there is a large radiation pressure 
force directed away from the X-ray source. The gas expands 
upwards (away from the disk) and sideways (along the disk) 
due to the combination of thermal pressure and radiation 
pressure. The gas is moved on distances Hs ^ H before the 
radiation force and illumination decreases so that the gas 
can both cool and fall back onto the disk. Note that the ma- 
terial cannot escape to infinity as it is gravitationally bound 
- the escape temperature from the inner accretion disk is 
T ~ 10" K (Begelman, McKee & Shields 1983). Note also 
that "local" winds do not occur for other than the magnetic 
fiare models (see, e.g., NKK; NK, §4.1), where the hydro- 
static balance assumption should be a good approximation. 
NKK approximated the effects of this local outflow by 
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introducing a "gravity parameter" A — J-g/Trad, which is 
the ratio of the vertical component of the gravitational force 
at the top of the disk, Tg = GM H p/R^ , to the radiation 
pressure force, J-rad = F-^atne/c. For a given F^, m and 
f corona, the gravity parameter is uniquely determined. How- 
ever, when the wind is present, hydrostatic balance does not 
apply as we discussed above. Therefore, NKK and NK al- 
lowed A to be a free parameter (even when _Fx is fixed) to 
account for the poorly constrained strength of the wind. The 
Thomson thickness of the skin, Ts, should always be lower 
if the wind is present compared with the values predicted 
by hydrostatic balance. By increasing the gravity parameter 
artificially, one also obtains smaller values of , hence qual- 
itatively one can account for the wind by making A larger 
than its true value. 

Here, instead of artificially increasing the gravity pa- 
rameter, A, we try to develop a simple approximation to a 
local wind that should place us in approximately the correct 
part of parameter space. In particular, we write 



where p is the gas pressure, and A is the characteristic scale 
on which the gas pressure is decreasing. If hydrostatic bal- 
ance assumption were appropriate, this scale would be the 
usual thermal pressure scale height (e.g., see A as defined 
in Nayakshin 2000). Here the geometry leads us to consider 
the pressure changing on scales of order the height of the 
flare, which is of order H (Nayakshin & Kazanas 2001b). 
The coefficient of proportionality can only be found in a full 
scale 2 or 3-D calculation of this problem. Simple estimates 
show it can depend on the ratio of the X-ray and disk fiuxes, 
Jx/fdisk, the Thomson thickness of the skin, disk accretion 
rate, etc. Therefore, we introduce a free dimensionless "wind 
parameter". A, to describe this uncertainty: A = AH. 

Simple estimates show that A 1 is appropriate to a 
thermally driven wind when the gas density decreases sig- 
nificantly at the distance ^ H of the fiare location. The case 
of small A corresponds to supersonic outflows, w S> Cs (cs is 
the sound speed), when A ~ Hcs/v. We will vary A below 
and study its influence on the reflected spectra. 

Summarizing, we use equation (^), together with A — 
AH, to replace the hydrostatic balance condition in the 
NKK code (eq. 19 in NKK). The lower boundary conditions 
- continuity in the gas density and T = Teff remain the 
same. The code then directly computes the density struc- 
ture of the X-ray irradiated disk including the effects of this 
outflow, so this is clearly an improvement on the previous 
description which used hydrostatic balance with the ad hoc 
gravity parameter, A. In terms of gas temperature profiles, 
the results are not appreciably different than those using A 
in that the nearly discontinuous transition between the hot 
and cold phases is still present. However, the advantage of 
the new approach is that we hope ^ automatically puts us 
in the roughly the correct parameter space for real magnetic 
flares. 



3 SIMULATED SPECTRA 

We simulate spectra over a range of expected parameters for 
the hard X-ray spectra seen from GBHC and AGN. Plainly 



the spectral index of the incident power law influences the 
results, so we test two values of the index; F = 1.6 and 
F = 1.8. The power law is assumed to extend out to 200 
keV, and then have an abrupt cutoff. This is a better model 
of a Compton scattered spectrum below ~ 80 keV than the 
very gradual rollover produced by an exponential cutoff. Ob- 
servations seem to indicate f corona ^ 0.8, so we fix this in 
all simulations. We assume an underlying disk accretion rate 
of m = 0.03 (this contrasts with the m = 10~^ assumed in 
paper 1) onto a super-massive black hole of 10* M©. The 
large black hole mass means that we avoid disk photons po- 
tentially contributing to the observed bandpass, and more 
importantly numerically we avoid the very high densities in 
GBHC disks for which the atomic data and codes are not 
reUable. XSTAR 2.0 (Kallman & Bautista 2001) has the 
highest density range of current photo-ionisation codes, but 
even this breaks down at ^ 10^* cm~^, whereas the GBHC 
disks can have densities of ~ 10^° cm"'' in the ionised skin, 
and 10^^ cm"'' in deeper layers. It is simply not possible to 
calculate the ion populations in X-ray illuminated GBHC 
disks with current codes. We use the AGN refiected spectra 
as the best current approximation to those from GBHC, but 
we caution that the GBHC ion populations can be different 
from those given here. 

To obtain fiat spectra, i.e. F < 2, requires that the 
power dissipated in heating the electrons in the flare is much 
larger than the seed photon luminosity available for Comp- 
ton scattering (e.g., Pietrini & Krolik 1995; Stern et al. 1995, 
Poutanen & Svensson 1996). Further, in order to produce re- 
flected spectra with little He-like Fe contribution, one needs 
the skin to be essentially completely ionized. Within a factor 
of few, both of these conditions require the X-ray flux inci- 
dent on the disk surface near flares footpoints to be much 
larger than the disk flux, i.e., F~^ 3> -Fdisk- We thus decided 
to test values of -Fx/-fdisk anywhere between ~ 10— 100, and 
we choose a grid of -Fx/^disk = 10, 30, 100. We also choose a 
grid for the wind parameter of A = 0.03, 0.1, 0.3 and 1. 

This results in 24 model spectra of reflection from a 
disk which extends down to the last stable orbit. We sim- 
ulate these spectra through the RXTE PGA response, as- 
suming a flux of ~ 10~* ergs cm~^ s^^ and a Iks exposure 
(i.e. comparable to a single PGA orbit of a bright GBHG). 
The resulting 3-20 keV spectra have systematic errors of 0.5 
per cent added in quadrature to the (negligible) statistical 
errors since it is sadly unlikely that any satellite will be cal- 
ibrated to better accuracy than this. These spectra are flt 
using the flxed density ionization models without relativis- 
tic smearing (pexriv: Magdziarz & Zdziarski 1995), with 
a separate narrow (intrinsic width flxed at 0.1 keV) Gaus- 
sian line of free normalization, and energy (except that this 
latter is constrained to be between 5.5 and 7.5 keV, as ex- 
pected for iron with some shifts allowed for Doppler and 
gravitational effects). All the flts were statistically accept- 
able ixl = 0.6 - 1.2 for 38 de grees of freedom) and gave 
a derived ionization consistent with mostly neutral mate- 
rial except for the case with Fx/-Fdisk = 100, A = 1. This 
extreme illumination with slow wind results in a large op- 
tical depth in the hot layer (> 1), so the observed reffected 
spectrum comes from the lower temperature material in the 
hot layer rather than the cool material underlying the skin. 
This gives rise to a clearly ionised reffected signature, so 
these data are not shown in the following plots. 
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Figure 1. Derived values of spectral index and reflected fraction for simulated RXTE dataiiles made from the illuminated disk models, 
fit with the pexriv constant density reflection model. The models are calculated for Fx/^diafc = 10 (large circles; upper panel), 30 
(crosses; middle panel), and 100 (squares; lower panel), each with wind parameter A = 0.03,0.1,0.3 and 1 (except for Fx/FiHgk = 100 
where A = 0.03,0.1 and 0.3). Stronger winds (small A) have smaller optical depth in the extremely ionized skin, and so have larger 
reflected fractions. The solid line in each panel is the observed index-reflection correlation of Zdziarski et al. (1999), while the dotted 
Une is that observed for GBHC (Gilfanov et al., 2000). 



Figure 1 shows the resulting derived photon spectral 
index, F, and solid angle, Q/2tt for each different Kc/Fdisk. 
For a given spectral index, the optical depth of the ionized 
layer is larger for higher Fx/fdisk, and for slower winds. The 
solid line shows the overlaid F — Q/2t: correlation derived 
by Zdziarski et al., (1999) for a sample of AGN and GBHC, 
while the dotted line shows that derived from GBHC by 
Gilfanov et al., (2000). Clearly, the simulated spectra can 
match the observed correlations only if the wind strength 
and Fx/Fdiak are varied in the "right way". Unfortunately 
values of these parameters depend crucially on the (currently 
unknown) physics of magnetic flares, e.g. on how the height 
of the flare might evolve as a function of mass accretion rate, 
so it is hard to see whether these parameters will indeed 
cooperate to produce the observed correlations. 

We get similar results (see Figure 2) in terms of derived 
spectral index and reflected fraction when using the reflec- 
tion code developed by Zycki et al., (1999). This code has 
the iron line intensity and energy tied self-consistently to the 
continuum reflection, and then both line and reflected con- 
tinuum are smeared by relativistic effects. Again the models 
fit the simulated data well, and imply a low ionization state 
of the reflecting material, but the difference now is that the 
spectral broadening of the line and reflected edge can be con- 
strained. Figure 2b shows the amount of smearing (in units 
of gravitational radii Rg = GM/ <? rather than Schwarzchild 
radii, Rs = 2Rg). In general a higher optical depth (slower 



wind and higher Fx/^disk) gives less smearing, but none of 
the spectra show the extremely smeared and skewed compo- 
nents which are associated with the last stable orbit around 
the black hole. 

Clearly the resolution from proportional counters leaves 
much to be desired. To investigate the smearing further we 
simulated the model spectra through the ASCA CIS re- 
sponse, assuming a mean 2-10 keV flux of ~ 10~* ergs cm~^ 
s~^, and a 20 ks exposure (as typical for the Cyg X-1 obser- 
vations: Ebisawa et al., 1996). We fit these spectra with the 
rclativistically smeared reprocessed spectral model of Zycki 
et al., (1999) over the 4-10 keV band (for comparison with 
the CIS Cyg X-1 data of Done & Zycki 1999). Again all 
the fits are adequate {xt = 0.75 — 1.25 for 79 degrees of 
freedom) with the reflected spectrum ionisation indicating 
mainly neutral material. Figure 3a shows the resulting de- 
rived spectral index and solid angle, showing that curvature 
in the data can be more easily matched in the narrower 
bandpass data by a somewhat steeper spectral index and 
correspondingly higher reflected fraction. Figure 3b shows 
the resulting inner disk radii. Again the reflected spectral 
features axe never as broad as expected if the reflector were 
a constant density disk extending down to the last stable 
orbit in a Schwarzchild geometry. 

This latter fact is somewhat surprizing given that the 
Thomson depth of the skin is quite low, Ts ~ 0.1, for the 
strongest wind models. Experimenting with fits, we found 



© 0000 RAS, MNRAS 000, 000-000 



X-ray irradiated disks 5 




Photon Index T Photon Index V 

Figure 2. (a) Derived values of spectral index and refleeted fraetion from the simulated XTE datafiles shown in Figure 1 but fit with 
constant density reflection models of Zycki et al., (1999). These include the sclf-consistently calculated iron emission lines for a given 
reflected continuum, and relativistically smear both line and eontiima together. The symbols are the same as for Figure 1. It is again 
clear that the optical depth of the extremely ionized skin is the major factor which determines the properties of the observed spectrum, 
but this is dependent on the (currently unknown) wind strength, (b) Derived values of the inner disk radius in units of Rg = GM/c^ 
obtained from the fits in (a) , plotted against the spectral index. While the RXTE resolution is fairly poor, the line is never as broad as 
expected if the reflector were a constant density disk extending down to the last stable orbit in a Schwarzchild geometry (= dRg). 
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Figure 3. As for Figure 2a and b, but with fits to simulated ASCA GIS datafiles. (b) The moderate resolution of the GIS detectors gives 

better constraints on the amount of relativistic smearing than those obtained from the RXTE simulations (figure 2b) . Plainly the ionized 
skin suppresses the line emission from the innermost regions, leading to derived disk radii which are always larger than 3 Schwarzchild 
radii (= 6Rg)- 



that the main reason that the inner disk radius is larger in 
the fits than in the actual models is that the broad redshifted 
wing of the Fe Ka lino gets confused with the incident power- 
law. Indeed, the inferred value of F differs from the actual 
one by as much as 0.1 for the strongest wind cases. 

4 REAL DATA 

Plainly these model spectra bear a close resemblance to the 
data in terms of the properties of the reflected spectrum. The 



next step is to fit such models to real data. As a first attempt 
at this we have taken the tabulated spectra used above and 
incorporated them into the XSPEC spectral fitting package 
(Arnaud 1996) as a local model. With this approach the 
only true free parameter is the normalisation of the reflected 
spectrum, but we have approximated the change made by 
small variations of spectral slope (AF < 0.1 by multiplying 
the models calculated for a given spectral index Fo by the 
factor E^~^'>. We have models tabulated for Fo = 1.6 and 
1.8, and these give reflected spectrum for F = 1.7 which 
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Figure 4. Data from the transient black hole Nova Musca taken on July 23 1991 when the source showed a classic low state spectrum. 
The data are well described by the complex ionization reflection models presented here. The dashed line shows the incident power law 
(which is absorbed at low energies by the interstellar column of 1.6 X lO'^^ cm~^) while the dotted line shows its reflection from a disk 



with a completely ionized skin (A 
spectrum. 



■ 0.3, Fx/Fdisii = 10, i = 63°). The lower solid line shows the contribution from the disk blackbody 



> 




channel energy (keV) 

Figure 5. ASCA GIS data for Cyg X-1 fit to the complex ionization reflection models presented here. The dashed line shows the incident 
power law while the dotted line shows its reflection from a disk with a completely ionized skin (A = 0.3, Fx/F^i^f^ = 100, i = 49°). The 
lower solid line shows the unsmeared, neutral reflection from the companion star/outer disk. 



are consistent to within 5% in the 2-20 keV band, so the 
total spectra for F = 1.7 derived from To = 1.6 and 1.8 are 
consistent to within 2 %. 

The first test is for broad band data from GINGA. A 

classic GBHC low state spectrum is that of Nova Muscae 
on July 23rd 1991. The constant density (single ionization 
parameter) reflection models give a best fit of = 15.3/24 
for a model consisting of a disk blackbody, together with a 
power law and its relativistic refiection component (inclinar 
tion fixed at 60°) with T = 1.72 ± 0.02, fi/27r = 0.24^°:^^, 
iiin = 50t^i?g, all absorbed by a (fixed) column of Nh = 
1.6 X 10^^ cm"^ (Zycki et al., 1998). We replace the power 
law and its relativistic reflection by the model spectra used 



above (incident fiux plus its refiection from the complex ion- 
ization structure of the disk) for To = 1.8 and step over all 
tabulated values of Fx/Fdiak and A. The best fit is shown in 
Figure 4, giving xl = 18.5/25 for F^/Fduk = 10,A = 0.3. 
Plainly low resolution, broad band data where the reflected 
spectrum appears to have much lower normalisation and 
smearing than expected for a disk around a black hole can 
be well fit by these model spectra which have the disk sub- 
tending a solid angle of unity and extending down to the 
last stable orbit. 

A higher resolution spectrum gives a better test of the 
relativistic smearing constraints. We use the ASCA GIS 6 
data from Cyg X-1 of Ebisawa et al., (1996), as analysed by 
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Done & Zycki (1999) in the 4-10 keV range. We step over 
inclinations from 56° —31° as the inclination of the Cyg X-1 
system is not well known, and include a narrow, neutral re- 
flected component (solar abundance, inclination of 60°, with 
line normalisation and energy tied to the reflected contin- 
uum) to account for the contribution to the reflected spec- 
trum from the companion star. The interstellar column is 
fixed at 6 X 10^^ cm"^ This gives a best fit of xl = 83.8/79 
for A = 0.3, F^/Fdisfc = 100, i = 49°, shown in Figure 5. 
There are many other solutions within Ax^ = 6.25 (90 per 
cent confidence for 3 parameters), although none of these 
allow the strongest wind with A = 0.03. The quality of the 
fit is plainly comparable to that of the constant ionization 
models of Done & Zycki (1999), who obtained xl = 81.4/79, 
n/2Tv = O.lOto ll and Ri^ = IStl^^g- Thus even high reso- 
lution data which appears to show much less reflection and 
(marginally) smearing from the accretion disk than expected 
can be fit with the models where the ionization instability 
gives rise to a completely ionized skin on the disk overlying 
more neutral material. 



5 CONCLUSIONS 

We use a modified version of the code of NKK and NK to 
calculate the vertical structure of an X-ray illuminated disk 
at a given radius, and then sum these spectra over all radii, 
with appropriate relativistic smearing, to get the full disk 
refiected spectra expected from magnetic flares. However, 
the refiected spectra are not uniquely determined as the local 
intense illumination expected from magnetic fiares can drive 
outflows from the illuminated region, lowering the optical 
depth of the ionized skin. We approximate the effects of this 
with a free parameter A (equation ^) . We simulate full disk 
reflected spectra for a grid of different Fx/Fdisk = 10,30 
and 100, with spectral index F = 1.6, 1.8 and with A = 
0.03,0.1,0.3 and 1. 

We simulate each total spectrum (illuminating power 
law plus full disk reflection) through the low resolution 
RXTE and moderate resolution ASCA GIS response and 
fit these with the constant density (single ionization state) 
pexriv models. The ionized skin reduces the derived Q./2ti 
for a given spectral index, and also reduces the amount of 
relativistic smearing observed as the ionized skin depth is at 
its maximum in the inner disk. 

The resemblance of the (low-to moderate wind veloc- 
ity) simulated spectra to observations of low-state GBHC 
and AGN is striking. The models have the disk sub-tending 
a solid angle of 2n and extending down to the last stable 
orbit, but the vertical ionization structure can mask the re- 
fiected signature especially from the inner disk. When such 
spectra are fit with the constant density refiection models 
then the derived solid angle and iron line width appear to 
imply that the disk subtends a small solid angle and has 
a large inner disk radius. These simulations of a full disk 
confirm the conclusions of DN and Ballantyne et al., 2001, 
drawn from local (single radius) models that the observa- 
tions of a small solid angle and line width do not necessarily 
rule out static magnetic fiares above an untruncated disk. 

We show for the first time that these magnetic flare 
models can fit real data. Both low resolution (GINGA data 
from the low/hard state of Nova Muscae, Zycki et al.. 



1998) and moderate resolution (ASCA GIS data from the 
low/hard state of Cyg X-1, Done & Zycki 1999) can be fit 
with these full disk models. The x^ for these fits is compa- 
rable to that obtained from the constant density reflection 
models which implied a truncated disk. Thus the data can 
be consistent with all the currently proposed X-ray pro- 
duction mechanisms - either simple (static) magnetic flares, 
or with outflowing plasma from a magnetic flare, or with a 
truncated disk/inner X-ray hot (advective) flow. 

To make progress we clearly need more sophisticated 
models of the reflected spectra. It is probable that the pri- 
mary reason that all these different models can fit the data 
is because of the current freedom in the values of model 
parameters (e.g., A for the wind from magnetic flares: the 
relativistic outflow velocity in the plasma ejection models: 
Beloborodov 1999; the transition radius in the advective flow 
models: Esin, McClintock & Narayan 1997). In particular, 
the magnetic flare models should be improved via an explicit 
calculation of the X-ray induced wind. This may limit the 
range of possible A for a given Fx/Fdisk, ruling some mod- 
els out. For the cold outer disk plus inner hot flow models, 
the next step is to include the presence of the skin on the 
cold disk - the physics of the thermal ionization instability 
applies to all configurations of the disk. Another modelling 
issue which needs to be addressed is the current inability of 
the atomic codes to handle the high densities required to 
properly model the GBHC disks. 

We also need better data, especially observations which 
extend over a wider bandpass. Including the higher energy 
spectra beyond 20 keV may start to break the degeneracy 
between the different model spectra (DN), while simultane- 
ous lower energy data can constrain the seed photon flux. 
Another possible way to distinguish between the models is 
to study their variability behaviour. The ionised skin of the 
disk can respond to changes in the illuminating flux on a 
timescale of the order of the Keplerian rotation time scale. 
This time scale is longer than the light crossing time at the 
same radius (Nayakshin & Kazanas 2001a). In other words, 
a highly variable light curve may produce reflected spectra 
that are different from the steady-state models (studied in 
this paper), which then could be used to differentiate be- 
tween the models. Finally, and perhaps most convincingly, 
Fe Ka line reverberation studies should present a variety 
of constraints on the models (Reynolds 2000, Revnivtsev et 
al. 1999). Recently, Nayakshin & Kazanas (2001b) showed 
that time-resolved Fe Ka line from a single magnetic flare 
is narrow but is moving across the ~ 4 — 8 keV band. This 
prediction is completely different from the instantaneous line 
profile expected in a central source or lamppost-like geom- 
etry (e.g.. Young & Reynolds 2000, Ruszkowski 2000), and 
one can hence hope that that will be the definitive test of 
the accretion flow geometry. 



6 ACKNOWLEDGEMENTS 

SN acknowledges support from NRG Research Associate- 
ship. The authors thank Demos Kazanas and Piotr Zycki 
for useful discussions. We are indebted to the anonymous 
referee for pointing out a very serious typographical error in 
the equation (ra). 



© 0000 RAS, MNRAS 000, 000-000 



8 C. Done & S. Nayakshin 



REFERENCES 

Arnaud K. A., 1996, Astronomical Data Analysis Software and 
Systems V, eds. Jacoby G. and Barnes J., pl7, ASP Conf. 
Series volume 101. 
Balbus S. A., Hawley J. F., 1991, ApJ, 376, 214 
Ballantyne D., Ross R., Fabian A. C, 2001, MNRAS, in press 
Begelman M. C, McKee C, Shields G., 1983, ApJ, 271, 70 
Beloborodov A. M. 1999, ApJ, 510, L123 

Chiang J., Reynolds C. S., Blaes O. M., Nowak M. A., Murray 
N., Madejski G., Marshall H. L., Magdziarz P., 2000, ApJ, 
528, 292 

Done C., Madejski G. M., Zycki P. T., 2000, ApJ, 536, 213 
Done C., Nayakshin S., 2001, ApJ, 546, 419 (DN) 
Done C., Zycki P. T., 1999, MNRAS, 305, 457 
Ebisawa K., Ueda Y., Inoue H., Tanaka Y., White, N. E., 1996, 
ApJ, 467, 419 

Esin A. A., McClintock J. E., Narayan R., 1997, ApJ, 489, 865 
Fabian A. C., Rees M. J., Stella L., & White N. E., 1989, MNRAS, 
238, 729 

Fabian A. C., Iwasawa K., Reynolds C. S., Young A., 2000, PASP, 

112, 1145 
Field G. B., 1965, ApJ, 142, 431 

Galeev A. A., Rosner R., Vaiana G. S., 1979, ApJ, 229, 318 
Gilfanov M., Churazov E., Revnivtsev M., 1999, A&A, 352, 182 
Gilfanov M., Churazov E., Revnivtsev M., 2 000. 5th CAS/MP G 



Workshop on High Energy Astrophysics (astro-ph/0002415) 
Haardt F., Maraschi L., Ghisellini G., 1994, ApJ, 432, 95 
Kallman T. R., Bautista M., 2001, ApJS, 133, 221 
Kallman T. R., White N. E., 1989, ApJ, 341, 955 
Ko Y.-K., Kallman T. R., 1994, ApJ, 431, 273 
Krolik J. H., McKee C. F., Tarter C. B., 1981, ApJ, 249, 422 
Lubinski P., Zdziarski A. A., 2001, MNRAS, 323, L37 
Magdziarz P., Zdziarski A. A., 1995, MNRAS, 273, 837 
Miller J. M., Fox D. W., Di Matteo T., Wijnands R., Belloni T., 

Pooley D., KouveUotou C, Lewin W. H. G., 2001, ApJ, 546, 

1055 

Narayan R., Yi I., 1995, ApJ, 444, 231 
Nayakshin S. 2000, ApJ, 534, 718. 

Nayakshin S., & Kallman T. R., 2001, ApJ, 546, 406 (NK) 
Nayakshin S., & Kazanas D., 2001a, ApJ, 553, L141 
Nayakshin S., & Kazanas D., 2001b, submitted to ApJ 
Nayakshin S., Kazanas D., & Kallman T. R., 2000, ApJ, 537, 833 
(NKK) 

Pietrini P., KroUk J. H., 1995, ApJ, 447, 526 
Poutanon J., Svensson R., 1996, ApJ, 470, 249 
Revnivtsev M., Gilfanov M., & Churazov E., 1999, A&A, 347, 
L23 

Reynolds C. S., 2000, ApJ, 533, 811 

Rozahska A., Czerny B., 1996, AcA, 46, 233 

Rozahska A., Czerny B., 2000, A&A, 360, 1170 

Ruszkowski M., 2000, MNRAS, 315, 1 

Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 

Stern B. E., Poutanen J., Svensson R,., Sikora M., Begelman M. 

C, 1995, ApJ, 449, 13 
Svensson R., Zdziarski A. A., 1994, ApJ, 436, 599 
Young A., Reynolds C. S., 2000, ApJ, 529, 101 
Young A., Fabian A. C, Ross R., Tanaka Y., 2001, MNRAS, in 

press 

Zdziarski A. A., Lubinski P., Smith D. A., 1999, MNRAS, 303, 
11 

Zycki P. T., Done C, Smith D. A., 1997, ApJ, 488, L113 
Zycki P. T., Done C, Smith D. A., 1998, ApJ, 496, L25 
Zycki P. T., Done C, Smith D. A., 1999, MNRAS, 305, 231 



